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^ IN^ _ . 

The generation of highly uniform magnetic fields is the primary goal of the design 
of magnets for NMR imaging in medical applications. The magnet described in 
this paper is a novel structure of permanently magnetized material that extends 
the field range to 0.5 Tesla, well above the current limit of commercially available 
permanent magnets, with a design approach that is consistent with size and 
weight well within practical limits. 

The magnet combines a conical structure of permanently magnetized material 
with ferromagnetic pole pieces. The geometry of the magnet is computed in Section 
2, that presents the principle of operation and the basic characteristics. The effect 
of the field perturbation due to the open magnet configuration is discussed in 
Section 3. 

2. CONICAL MAGNET 

Fig. 2.1 shows the basic conical structure of magnetized material contained between 
two coaxial cones of half angles e r e Q . The material is assumed to be magnetized 
with a uniform remanence J Q parallel to the axis Z of the cones. The medium 
surrounding the external cone is assumed to be non-magnetic, and the medium 
inside the internal cone is assumed to be a ferromagnetic material with infinite 
magnetic permeability. 

Assume a system of spherical coordinates n, e, 0 where n is the distance of a 
point from the center O of the cones, e is the angle between n and the axis z of 
the cones, and 0 is the angle between the plane formed by a point P and the z 
axis and an arbitrary plane that contains the Z axis. 
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Report [1]. Assume the limit of infinite radial dimension of the cones and consider 
the special case of a field configuration independent of the angular coordinate 0 . 
Let 



x = cos h 



(2.1) 
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Fig. 2.1 Conical Structure. 
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Fig 2.3 Equipotential lines. 



The magnetostatic potential 0 (n, b) generated by remanence J 0 within the 
magnetized material is given by the particular solution of the Laplace's equation 
that satisfies the boundary condition 



0(n, e) = 0 



(2.2) 



One has 



(p(p,e)=-^(i-xj)xo P 



X. 2 



(1+X,)(1-X) 



(1-X,)(1+X) 



(2.3) 



for e, < e < e 0 and 6 - e 0 < e < 6 - e, i 0 is the magnetic permeability of a 
vacum and 



x^cose,,, x 1 = cosej 
In the region e fl < e < 8 - e fl the magnetostatic potential is 



d>(p,e) = -A(i_ x ^xp 

M-o 



1 X ° I ^ (^H 1 ***) 



Xl X2 (l-Xl)(l-Xo) 



(2-4) 



(2.5) 



Eq. (2.5) is the potential of a uniform magnetic field of intensity H 0 oriented parallel 
to the z axis and of magnitude 



1 0 H Q - K J 0 



where 



x, 2 (1 -*,)(! +x 0 ) 



(2-6) 
(2.7) 



The equipotential surfaces are planes perpendicular to z in the region e 0 < e < 8- 
e 0 and within the magnetized material the equipotential surfaces are given by the 
equation 



P = 



2 (l-xi)(l + x) xi 



(2.8) 



where 

The field inside the cone is not uniform and the values of H and B at the surface 
of the cone e = £ 0 are given by the vector diagram of Fig. 2.2 

An example of the equipotential lines given by Eq. (2.8) is shown in Fig. 2.3 for 
the cone angles c l = 20°, e 0 = 45°. In this example, the uniform field in the 
region e Q < e < 6 - e Q corresponds to a value of K as given by Eq. (2.7) 

K = 0.426 (2.10) 

Hence if one assumes a value 

J 0 - 1.3 Tesla (2.11) 
the cone generates a uniform field of intensity H 0 

i 0 H 0 - .55 Tesla (2.12) 

The basic schematic of Fig.2.1 can provide the computation of the field in a 
structure where the magnetized material is eliminated in the region 

-Z 0 < Z < + Z 0 (2.13) 

The new structure is shown in Fig. 2.4 where the two areas of the surfaces that 
limit region (2.13) are circles of radius 

r 0 =z 0 tan£ 0 (2.14) 

The two circular areas are assumed to be surfaces of an infinite magnetic 
permeability material. Furthermore, the equipotential surface within the magnetized 
material that start from the circle of radius r o is also assumed to be transformed 
into a surface of infinite magnetic permeability material. As a consequence, the 
field configuration inside the magnetized material is still given by Eq. (2.3) and the 
structure of Fig. 2.4 generates the uniform field defined by Eqs. (2.6),(2.7) 
everywhere outside the magnetized material, including the gap between the 
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Fig. 2.5 Basic schematic of the magnet. 




2.6 Computation of the remanence of the transition region (ABC). 



ferromagnetic plates of radius r 0 . An additional source, other than the magnetized 
material in Fig. 2.4, is necessary to maintain the block of ferromagnetic material 
in the structure of Fig 2.4 at the potential 




(2.15) 



The design of the complete structure depends upon the external termination of 
the conical body of magnetized material. This can be accomplished by the transition 
from the conical structure to a cylindrical structure shown in the schematic of Rg. 
2.5. The cylindrical structure is contained between the planes Z = Z v z = z 2 and 
its thickness Z 2 - Z 1 is given by 



The cylindrical structure is magnetized with the uniform remanence J 0 parallel to 
the axis Z. The structure cannot interface directly with the cone and a region (ABC) 
of magnetic material of remanence different from J has to be inserted between 
the cone and the cylindrical structure, as shown in Fig. 2.5. The computation of 
the geometry of region (ABC), and its magnetization is complicated by the cylindrical 
configuration of the magnet. However, even in the approximation of a two- 
dimensional geometry independent of the coordinate perpendicular to the plane 
of the section of Rg. 2.5, the transition between the two median would require a 
uniform remanence oriented along the z axis of magnitude [2] 



as shown in the vector diagram of Rg. 2.6, where H, B indicate the intensity and 
the induction in the cone, and H 1 , B 4 indicate the intensity and the induction in the 
cylindrical structure. Because the intensity in the transition region coincides with 
the intensity in the cylindrical region the interface (AB) between the cone and 
region (ABC) is a line w of the points P of intersection of the equipotential lines 
O in the cone with the same equipotential lines within (ABC). 



Z2 + Z1 
Z2 



(2.16) 



J =2J 0 



(2.17) 



By virtue of Eqs (2.8) and (2.9), the equation of the equipotential lines inside the 
cone can be written in the form 



xa rn (lH-xi)(l-xo) t t xo 

p 2 (l-xi)(lH-xo) xi 

P~ x /n (l + xQ(^x) x (2.18) 

2^(l-xi)(l+x) + xi 

where ft is the distance of the point of potential 0 on the cone of angular width h 0 
from 0. In the coordinates it, e, the equations of the equipotential lines in the 
region (ABC) are 

Z Z2-Z (D 

p = -, = - (2.19) 

X Z2-Z1 <pi 

Thus, the system of Eqs. (2.18) and (2.19) provides the equation of line w 

xo £n (l+*0(l-xo) | t xo 
l_p 2 (l-x.)(l + xo) x. (220) 

x z 

2 (l-xi)(l + x) xi 

The potential 0 of the i = <*> component of the magnetic structure of Fig. 2.5 is 
generated with the additional planar structure contained between the planes Z = 
Z Z = z where Z 0 can be chosen somewhat arbitrarily, and the thickness 

3, 4 3 

z 4 -z 3 is given by 

Z *~ Z * =K (2.21) 

Again, this externa! planar structure is uniformly magnetized with the remanence 
J 0 parallel to the z axis. The surface z = z 4 is the interface between the magnetic 
structure and the external yoke of the magnet that closes the flux of the magnetic 
induction. 
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3.1 Equipotential lines with no transition structure between cone and cylinder. 



r /r=2. THETA =45. THETA 1 =25, long strap 
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Fig. 3.2 Equipotential lines in the transition region with the insertion of the 
three straps. 




Fig. 3.3 Held on a reference sphere with center at 0 and diameter 2z 0 ,(a) 
conical magnet, (b) conventional magnet. 




Fig. 3.4 Equipotentlal lines with a single strap and K*= 
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Fig. 3.5 Field over the reference sphere for increasing valuve of 



3.DESIGNC0NS!DER/^NS. 



The exact solution presented in Section 2 is affected by two main constraints in 
its implementation in a practical design. In the first place, a finite external radius 
has to be chosen and its value has to be large enough to minimize the perturbation 
of the field uniformity in the gap between the pole pieces. Second, the design 
presented in this paper is aimed at the development of a magnetic structure 
capable of generating the largest possible field compatible with practical 
considerations. This means that the highest commercial V available remanence J 
would be chosen for the magnetic material. If a value J 0 ~ 1.3 T is chosen, as 
assumed in Section 2. no permanent magnet of remanence 2 J 0 is available for 
the implementation of the exact designs of the transition between the cone and 
the cylindrical sections. A remanence 2 J 0 is well above the theoretical limit of 
what can be achieved in the development of rare earth magnetic materials. As a 
consequence an approximate solution may be adopted taking advantage of the 
fact that the transition region is well outside the region of interest is the gap 
between the pole pieces. The approximate solution adopted in this design 
eliminates the transition region (ABC) and assume a uniform distribution of J 0 
throughout the region of the cone and the cylindrical region as shown in Fig. 3.1. 
The effect of the field distortion is illustrated in the example of Fig. 3.1 which 
shows the result of the numerical calculation of the equipotential lines. A partial 
correction of the field distortion in Fig. 3.1 is possible with the insertion of 
magnetized straps within the transition region. The straps are relatively thin plates 
of high magnetic permeability material deployed along the theoretical equipotential 
surfaces within the transition region. Schematically, the insertion of the straps is 
shown in Fig 3.2 where the straps are extended over the entire radial dimension 
of the cylindrical structure. 

Ideally, the region of interest should extend to the full region of the gap between 
the pole pieces. Consider the value of the field intensity over a sphere whose 
center coincides with the center of the gap and whose diameter is equal to the 
dimension 2z 0 of the gap. An example of the numerical calculation of the field 
intensity over circles on the sphere in planes perpendicular to the axis z js shown 
in Fig. 3.3. One observes that the field at the points z = ± Z 0 is lower than the 
field on the equational plane Z = 0. 



0 

K relative to nominal (%) 



Fig. 



3.6 Field homogeneity and field intensity versus K' 




;. 3.7 View of the magnet with the external yoke. 



This is contrary to what is normally found between the pole pieces of a conventional 
magnet, where the filed increases from the center along the axis and decreases 
with the radial distance in the plane perpendicular to the axis as shown in Fig.3.3. 
Thus, in the implementation of the conical magnet, the field generated by the 
cylindrical region over compensates the field generated within the central gap. As 
a consequence, the external planar structure contained between the planes z 3 < 
z < z 4 can be designed for a value K' > K and the value of K' can be chosen to 
minimize the variation of the field intensity in the region of interest. 

The results of the selection of K" are shown In the magnet configuration of Fig 3.4 
where a single strap is used for the partical correction of the field distortion in the 
transition between cone and cylindrical structure. The equipotentiai lines in Fig 
3.4 correspond to 

K' =1.03K ( 31 > 

The plotting of the intensity over the reference sphere of diameter approximately 
equal to 2z 0 is presented in Fig. 3.5 which shows that for K' largerthan the value 
(3.1) the field intensity on the equatorial plane becomes larger than the intensity 
on the z axis. At the optimum condition the field homogeneity on the reference 
sphere is better than ~ 1 %. The homogeneity changes rapidly with differentvalues 
of K' as shown in Fig. 3.6. As expected, the field intensity within the sphere 
increases as K' increases. As shown in Fig. 3.5 the field at the center of the 
sphere is 

j o h q ~ .0.47 Tesla (3-2) 

for J 0 ~ 1.38 Tesla. The homogeneity at condition (3.1) is a remarkably small 
value as a design parameter of the open magnet presented in this paper, and one 
can expect the tolerances of fabrication and material properties to have a larger 
affect on a practical implementation of the design approach. Hence, a conventional 
shimming technique can be used to improve the homogeneity to a value of the 
order of 50 ppm or better consistent with imaging requirement. 



The schematic of Fig. 3.4 corresponds to the basic parameters 



K = 0.51, e 0 = 45°, e, = 16° 
r-.-l.75r, 

^COMPE NSATION O F FIEL D PISTOgrjON 

Assume that the region of interest is the sphere of diameter 2z 0 and center located 
at the center 0 of the magnet, where the field uniformity must satisfy the imaging 
requirements. A distortion of the uniform field is caused by the finite dimensions 
of a practical magnet and by the fabrication tolerances of the geometry and 
distribution of the magnetization. As shown in Section 3, a partial computation of 
the effect of finite dimensions can be achieved in the design phase by modifying 
the value of K of the external magnetic structure. The distortion due to the 
tolerances cannot be compensated in the design phase and requires the use of a 
shimming technique based on the measurement of the field generated by the 
assembled magnet The shimming technique requires the insertion of a magnetic 
structure whose configuration is dictated by a logic based on the field 
measurements. 

Assume that the shimming is performed by means of a distribution of magnetic 
dipoles on the conical surface (P Q B) of the schematic of Fig. 2.5 as shown in Fig. 
4.1. The dipole distribution is defined at each point P by a dipole moment density 
m per unit surface of components m r> m z as indicated in Fig. 4.1. Assume that 
m , m z are independent of the angular coordinate ^ . The dipole distribution m z 
on the element of surface 2itrdpof a ring of radius r coaxial with the axis Z 
generates a field on the sphere of radius p 0 shown in Fig. 4.2, whose potential is 
shown in Fig 4.3 for six equally spaced positions on the element of surface line 
(P 0 B) The z component of the intensity corresponding to the six curves of Fig. 4.3 
is shown in Fig. 4.4. The effect of the orientation of the dipole moment density is 
Illustrated in Fig. 4.5 where the z component of the intensity is computed for the 
component m r of the dipole moment density. The fields generated by m r , and rr\ 
exhibit a similar behavior. In both cases the field on the reference sphere decreases 
rapidly with the distance from point P Q . The significant difference between the two 
cases is the larger intensity generated by m r located in proximity of P 0 where m r 
is essentially perpendicular to the surface of the ferromagnetic pole piece. 



(3.3) 



Assume that the chosen dipole distribution is limited to the component m r . The z 
component of the intensity on the reference sphere of radius p 0 can be written in 



the form 

Pi 



Po (4.D 



Where y (p,6) is a function of the geometry of the magnetic structure and the 
distance between a point of the reference sphere of coordinates e vand a point 
of the conical surface 8 = 6. with the same angular coordinate V . The complexity 
of the geometry of the magnetic structure makes it impossible to achieve an 
exact analytical solution of Eq. (4.1) and numerical calculations have to be used 
to map function y as a function of 8 for each value of p . The maximum range of 
interval in Eq. (4.1) extends to the length P 0 B of F.g 2.5. With the assumption 
that the dipole distribution in confined to n+1 rings equally spaced alone the 
radial coordinate p at a distance 

g p = PaZPo (4.2) 

n 

Eq. 4.1 transforms to 

M- 0 Hz(e) = E tfMMfO (4.3) 



where 

P (Ph) = ™ r W Sp ( 4 - 4 > 



p 0>Pn are the radial coordinates of points P 0> B respectively. The field on the 
reference sphere can be expanded in spherical harmonicics 



^Hz(e) = I> k PK(cose) 

k=0 



(4.5) 




Fig. 4.1 Compensating dipole 




4.2 Schematic for the computation of the dipole field 




Fig. 4.3 Potential on the reference sphere generated by a distribution of dipoles 
oriented along the axis z. 



28 




Fig. 4.4 z component of the field generated by the z component of the dipole 
distribution 




Fig. 4.5 z component of the field generated by the r component of the dipole 
distribution 



Where P k is the Legendre's polynomial of the first kind defined as 



1 2 k k!dx^ ; (46) 

X = COS0 



if the field satisfies the symmetry condition 

H z (*-0) = H (9) (4.7) 

the expansion (4.5) is limited to the even values of k. The coefficients b k of the 
expansion are given by 

b k = ^EC, h P(pJ (4.8) 

where 

+1 

C k , h= = / ^(p hf x)P k (x)dx (4 9) 

-1 

Identically, the distortion § Hz of the field generated by the magnet on the reference 
sphere can be expanded in the spherical harmonicics 

oc 

Mo 8H(x) = $>P,,(x) (4.10) 



k=0 



where the coeffecients a k are 



a K = ^(2 k +l)J[6h z (x)]P k (x)dx (4 . 13L) 
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Fig. 4.6 Intensity on the reference sphere for K' = K and K' = 1 . 03 




Fig. 4.7 Partial compensation of the field distortion with two dipole rings 




Fig. 4.8 Schematic of a two ring dipole structure 



and the field distortion is eliminated if coefficients b k satisfy the system of 
equations 



(k = 0,1,2,...) 



(4.12) 



i.e. the distribution of dipoles on the conical surface G = Go that cancels the field 
distortion is given by the solution of the system of equations 



In a practical situation the compensation of the field distortion can be confined to 
a limited number of the lowest order harmonies, which as a consequence results 
in the same limited number of rings of the compensation structure of dipoles. 

The compensation with the dipoie distribution on the conical surface 6 = Bo can 
be illustrated in the simple example of the magnetic structure introduced in Section 
3 and depicted in Fig. 3.4. Consider the case of the dual distribution of magnetic 
material where the value of K' of the external layer in the region in larger than the 
value in the rest of the magnetic structure. Fig. 4.6 show the plotting of the 
intensity over the reference sphere for K'= Kand for the optimum value (3.1) of K* 
which minimizes the field distortion within the sphere. A common zero level of the 
field in Fig 4.6 is arbitrarily chosen for the two values of K\ In the particular case 
of K 4 = 1.03K the zero level corresponds to an average value of approximately 

lesia aime cenitn ui me syndic as ohovvm in • igt.u uio now auamo 
same value of 0.47 Tesla at both ends of the angular integral e = Oand 9 =90° 
. In the structure selected for the computation of Fig. 4.5, assume that the dipoie 
distribution is limited to the first two rings, p = P 0 » P = Pi i.e. 




n 



+1 



(k =0,1,2,...) 



p(Ph) = 0 
h>l 



(4.14) 



Select the values of p(p 0 )and p^that satisfy the equation 




(4.15) 



The total field generated at =0 by the two rings is 



m>S±H 2 (0) = p(Ri)¥(M>) + p(p 0 )^(po,0) 



(4.16) 



Thus the combination of the fields generated by the two rings over the reference 
sphere is the sum of a term independent of 0 given by Eq. (4.16) and a second 
term 6^(8) that is zero at 9 = 0 and 9 = 90° and has a maximum as shown 
in Fig. 4.7. 

The same figure shows the field generated by the two rings 



In the first approximation the correction of the distortion is accomplished by 
selecting the value of the dipole moment distribution that makes the maximum of 
\x$ 2 Hz ec l ual t0 ^ e maximum of value computed in Fig 4.5 for K 1 = 

1.03 K. As apparent from Fig. 4.5 the residua! distortion is of the order of 
10 3 ppm. 

The term generated by the rings does not contribute to the field distortion and its 
effect is a small change of the field of the magnet. A schematic of the rings 
configuration is presented in Fig. 4.8 
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